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Abstract 

Starting from the well-known field theory for directed percolation, we describe an 
evolving population, near extinction, in an environment with its own nontrivial 
spatio-temporal dynamics. Here, we consider the special case where the environment 
follows a simple relaxational (Model A) dynamics. Two new operators emerge, with 
upper critical dimension of four, which couple the two theories in a nontrivial way. 
While the Wilson-Fisher fixed point remains completely unaffected, a mismatch of 
time scales destabilizes the usual DP fixed point, suggesting a crossover to a first 
order transition from the active (surviving) to the inactive (extinct) state. 
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Introduction. The understanding of ecological catastrophies, where one or sev- 
eral species become extinct, is a major challenge in various areas of science. 
From a physics perspective, extinction events are frequently associated with 
a continuous phase transition from an active to an inactive (absorbing) state 
from which the population cannot recover. As the transition is approached, 
the dynamics is characterized by large fluctuations, whose characteristic length 
scale diverges. By virtue of this diverging scale, the large- distance long-time 
behavior of systems near continuous phase transitions is universal, i.e., inde- 
pendent of microscopic detail, and can be captured successfully by minimal 
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models. Renormalization group (RG) approaches have been extremely suc- 
cessful in exploring the physics of such models near criticality. In this brief 
note, we explore the effects of a fluctuating environment, reflecting, e.g., a food 
supply with its own nontrivial dynamics, on the universal properties of this 
active-inactive state transition. We first summarize the field theory for this 
transition in a uniform environment. Next, we introduce the dynamics of the 
environment and discuss its effects on the evolving population. We conclude 
with a brief summary and some open questions. 

The model. The minimal model describing the extinction of a single species in 
a uniform environment is well known [1,2,3]. It finds its applications in a broad 
range of problems, ranging from population dynamics to catalysis, forest fires, 
and directed percolation (DP). It is most easily expressed in the language 
of "chemical" reactions, describing the birth (^4 — > 2A), death (^4 — > 0), 
and overcrowding {A-\- A ^ A) of a population, with rates a, /x, and g/2, 
respectively. The individuals are free to diffuse in d dimensions. Once the last 
individual has died, the species cannot recover. Denoting the space- and time- 
dependent density of the population by s(a;, t), the associated field theory is 
well known [1,2] and can be written in terms of a dynamic functional, with a 
response field s{x,t): 

Jdp{s, s}^J d^x J dtXs{X-'dt + - V') + ^ {gs - gs)}s (1) 



Here, t = /j, — a is the net death rate which controls the distance from the 
(mean-field) critical point, and A sets the time scale. The stochastic character 
of the dynamics, along with the absorbing state condition, is refiected in the 
noise vertex, gs'^s. Invariance under rapidity reversal implies g = g. For now, 
we retain different nonlinearities here, to allow for a breaking of this symmetry 
in the full model. The effects of the environment are encoded in the constant 
rates r, g, and the time scale A. It is natural to study the consequences of 
fiuctuating rates, refiecting an environment that has its own nontrivial spatio- 
temporal structure. Previously, such effects as a quenched random t{x) [4] or 
a coupling to a diffusive mode [5,6,7] were investigated. Here, we study an 
environment with non-conserved relaxational dynamics, described by Model 
A [8,9]. Denoting the Model A local order parameter as ip, modelling, e.g., a 
nutrient supply for the population, the corresponding field theory is given by 



Ja{4>, 4>} = j d'^x J dt^{ip ^-'dtip + (r - v') <^ + |<^' 



- (2) 



where ip is, again, the corresponding response field, and r is the critical param- 
eter. Noting that both field theories have an upper critical dimension dc = 4, 
we seek novel couplings, involving both Model A and DP fields, which might 
destabilize the famihar DP [2] or Wilson-Fisher [10] fixed points. If /i is an 
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external momentum scale, naive dimensional analysis results in r ~ r ~ /i^. 
Assuming that the net death rate depends on the environment (e.g., via the 
availability of nutrients) and that the presence of the population affects its 
environment, it is natural to write an expansion t — ro + + wip^ + ... and 
r — To -\- vs -\- .... Dimensional analysis shows that k ~ v ~ 

and w ~ iJi,^~'^. The coupling k is relevant in = 4 and must be tuned to zero to 
access the multi-critical point r = r = k, = 0. Alternately, if we demand that 
the up-down symmetry of the Ising model remain valid, k can be set to zero 
for physical reasons. Once zero, it is not generated under the RG. Collecting, 
our model is given by J{0, 0, s, s} = Jdp{s, s} Ja{0, 0} + Jint{4>, <t>, s, s} 
where 

Jint{4>,4>,s,s} = j (fx j dt ^^v(j)(l)s + ^Awss^^l (3) 



Expectation values of the four fields are given as functional integrals with 
weight exp (—J). 

RG results. All five nonlinear couplings - u, g, g, v, and w - are marginal in the 
upper critical dimension dc — 4:. The new vertex v violates the invariance under 
rapidity reversal but respects the absorbing state condition. We use renormal- 
ized perturbation theory, in e = dc — d, combined with minimal subtraction. 
Due to the absence of a bare correlator in the DP field theory, no corrections 
to the Model A correlation and response functions are generated. As a result, 
the Model A fixed point remains at the Wilson-Fisher value, and all Model A 
exponents retain their familiar values. There are, however, corrections to the 
DP couplings, and nontrivial fiow equations for the new mixed couplings. In 
one-loop order, we find numerous infrared unstable fixed points, and one non- 
trivial stable one, characterized by the renormalized values u = 2e/3 + O(e^), 
gg = 4e/3 + 0(e2), vg ^ + 0{e^), and w = 13 (A + 7) / (12A). A geometric 
factor Gg = 2r (1 + e/2) (An)"'^^^ occurs in each loop integral and has been 
absorbed in the definition of the couplings. At first glance, this looks encour- 
aging; in particular, vg = restores the rapidity reversal symmetry. However, 
at this fixed point, the ratio p = 7/A of the two time scales fiows towards 
p — 0. This implies that the dynamics of Model A freezes on the time scale 
of the DP fields. The limit p — > turns out to be singular, so that the theory 
needs to be be reanalyzed. At the most naive level, this involves adding the 
time-delocalized vertex, w J d'^x J dt J dt' s{x,t)s{x,t)s{x,t')s{x,t'), to the DP 
field theory. Once included, however, all fixed points become unstable [3,4], 
indicating that the active-inactive state transition might turn first order. At a 
more rigorous level, one should integrate out the Model A fields and then take 
the static hmit. The resulting field theory will be analyzed elsewhere [11]. 

Conclusions. In summary, we have investigated the effect of a fiuctuating 
environment on the evolution of a population near the brink of extinction. 
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These fluctuations may be due to a food source, for example, which has its own 
nontrivial spatio-temporal dynamics. Here, we consider the interactions of the 
well-known DP field theory with an auxiliary field obeying simple relaxational 
(Model A) dynamics. Investigating the combined field theory to one-loop order 
in e = 4 — d, we find that the RG fiow implies a freezing of the Model A 
dynamics on the times scales of the DP fields. Reanalyzing the theory in this 
(singular) limit, all fixed points are now infrared unstable, indicating a possible 
first order transition. Two open questions remain, namely, first, to quantify 
the crossover to the time-delocalized theory, and second, to compare these 
field-theoretic findings near dc — 4: to Monte Carlo simulations in two and 
three dimensions. 
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